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ABSTRACT 

We present a detailed study of a set of gas-dynamical simulations of galaxy groups and 
clusters in a flat, ACDM model with O m = 0.3, aimed at exploring the effect of non- 
gravitational heating on the observable properties of the intracluster medium (ICM) . 
We use GASOLINE, a version of the code PKDGRAV that includes an SPH description of 
hydrodynamics to simulate the formation of four cosmic halos with virial temperatures 
in the range 0.5;$ T^S 8 keV. These simulations resolve the structure and properties 
of the intra-cluster medium (ICM) down to a small fraction of the virial radius, i? v ir- 
At our resolution A-ray luminosities, (Lx), of runs with gravitational heating only 
are in good agreement, over almost two orders of magnitude in mass, with analytical 
predictions, that assume a universal profile for CDM halos. 

For each simulated structure, non-gravitational heating of the ICM is implemented 
in two different ways: (1) by imposing a minimum entropy floor, Sfi, at a given redshift, 
that we take in the range 1^5 z ^5; (2) by gradually heating gas within collapsed 
regions, proportionally to the supernova rate expected from semi-analytical modeling 
of galaxy formation in halos having mass equal to that of the simulated systems. 

Our main results are the following, (a) An extra heating energy Ehlt I keV per 
gas particle within R V i r at z = is required to reproduce the observed Lx~T relation, 
independent of whether it is provided in an impulsive way to create an entropy floor 
Sfi = 50-100 keV cm 2 , or is modulated in redshift according to the star formation 
rate; our SN feedback recipe provides at most Eh — 1/3 keV/part and, therefore, its 
effect on the Lx~T relation is too small to account for the observed Lx~T relation. 
(b) The required heating implies, in small groups with T ~ 0.5 keV, a baryon fraction 
as low as is 40% of the cosmic value at i? v j r /2; this fraction increases to about 80% 
for a T ~ 3 keV cluster, (c) Temperature profiles are almost scale free across the 
whole explored mass range, with T decreasing by a factor of three at the virial radius. 
(d) The mass-temperature relation is almost unaffected by non-gravitational heating 
and follows quite closely the M cx T 3 / 2 scaling; however, when compared with data 
on the MsoQ-T ew relation, it has a ~ 40% higher normalization. This discrepancy is 
independent of the heating scheme adopted. The inclusion of cooling in a run of a 
small group steepens the central profile of the potential well while removing gas from 
the diffuse phase. This has the effects of increasing T ew by ~ 30%, possibly reconciling 
the simulated and the observed MsoQ—T ew relations, and of decreasing Lx by ~ 40%. 
However, in spite of the inclusion of SN feedback energy, almost 40% of the gas drops 
out from the hot diffuse phase, in excess of current observational estimates of the 
amount of cold baryons in galaxy systems. 

Likely, only a combination of different heating sources (SNe and AGNs) and cool- 
ing will be able to reproduce both the Lx~T ew and M^-T ew relations, as observed 
in groups and clusters, while balancing the cooling runaway. 
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1 INTRODUCTION 

The X-ray emission from clusters of galaxies offers a unique 
means for studying the physics of cosmic baryons and their 
connection with the processes of galaxy formation and evo- 
lution. The first attempt to model the ICM in the framework 
of the hierarchical clustering scenario assumed its thermo- 
dynamical properties to be entirely determined by gravita- 
tional processes, such as adiabatic compression during col- 
lapse and shock heating by supersonic gas accretion (Kaiser 
1986). Since gravity does not have characteristic scales, for 
an Einstein-de-Sitter cosmology this model predicts rich 
massive clusters and poor groups to appear as scaled ver- 
sions of each other. Under the assumptions of emissivity 
dominated by free-free bremsstrahlung and of hydrostatic 
equilibrium of the gas, this model predicts Lx oc T 2 (l+z) 3 / 2 
for the shape and evolution of the relation between X-ray 
luminosity and gas temperature. Furthermore, if we define 
the gas entropy as S — T/n 2 J A (n e : electron number den- 
sity; e.g. Eke et al. 1998), then the self-similar ICM has 
5ocT(l + 2 )- 2 . 

This simple model was quickly recognized to fail at ac- 
counting for several observational facts: (a) the Lx~T re- 
lation for nearby clusters is steeper than predicted, with 
L x oc T~ 3 for T> 2 keV clusters (e.g., David et al. 1993, 
White, Jones & Forman 1997, Allen & Fabian 1998, Marke- 
vitch 1998, Arnaud & Evrard 1999), with a possible further 
steepening at the group scale, T^z 1 keV (Ponman et al. 
1996, Helsdon & Ponman 2000); (b) no evidence of evolu- 
tion for its amplitude has been detected out to z<) 1 (e.g., 
Mushotzky & Scharf 1997, Donahue et al. 1999, Fairley et al. 
2000, Delia Ceca et al. 2000, Borgani et al. 2001b, Holden 
et al. 2002); (c) the gas density profile in central regions 
of cooler groups is relatively softer than in clusters and, 
correspondingly, the entropy is higher than predicted by 
self-similar scaling (e.g., Ponman, Cannon & Navarro 1999, 
Lloyd-Davis et al. 2000, Finoguenov et al. 2002); (d) colder 
clusters contain a relatively smaller amount of gas than very 
hot ones (e.g. Neumann & Arnaud 2001). 

A common interpretation for such observational facts 
requires non-gravitational energy input, which should have 
taken place during the past history of the ICM. This extra 
heating would place the gas on a higher adiabat and sustain 
the gas during the gravitational collapse of the cluster dark 
matter (DM) halo, preventing it from reaching high density 
in central regions and suppressing the X-ray emission. For 
a fixed specific amount of heating energy per gas particle, 
the effect is larger for the smaller systems, while being neg- 
ligible for more massive, hotter systems. Kaiser (1991) sug- 
gested for the first time that gas pre-heating and subsequent 
adiabatic collapse turns into a differential steepening of the 
Lx~T relation at the cluster scales. Evrard & Henry (1991), 
Navarro, Frenk & White (1995) and Bower (1997) assumed 
the gas in the cluster core to have a minimum entropy level, 
which was established at some pre-heating epoch. 

Cavaliere, Menci & Tozzi (1998, 1999) were able to pre- 
dict the correct Lx~T slope, both at the cluster and at the 
group scale, by assuming a pre-heated gas at a tempera- 
ture T ~ 0.5 keV, which is shocked with different strengths 
when falling into groups and clusters. A pre-heating pro- 
ducing an isentropic gas distribution has been assumed by 
Balogh, Babul & Patton (1999) and Tozzi & Norman (2001, 



TN01 hereafter). In particular, TN01 worked out a series of 
predictions of observables properties of the ICM with isen- 
tropic pre-heating, including the amount of energy feedback 
required to produce the correct Lx~T relation and entropy 
threshold at the group scale. Brighenti & Mathews (2001) 
discussed the effect of heating the gas either when it has 
still to collapse within clusters and groups (external heat- 
ing) or when it is already accreted (internal heating) . They 
concluded that internal heating is more efficient at account- 
ing at the same time for both the slope of the Lx~T relation 
and the excess entropy in poor clusters. 

While most of such analysis agree that an extra heating 
energy of ~ 0.5-1 keV per gas particle is required, yet no 
general consensus has been reached about its astrophysical 
origin. Supernovae (SNe) have been advocated by several au- 
thors as a possible source for pre-heating (e.g., Wu, Fabian 
& Nulsen 1998) and metal enriching (e.g. Loewenstein & 
Mushotzky 1996, Renzini 1997, Finoguenov & Ponman 1999, 
Pipino et al. 2002) the ICM. Although some analyses showed 
that SNe can actually provide an adequate energy budget 
(Menci & Cavaliere 2000, Lloyd-Davis et al. 2000), other 
authors claimed that this requires a very high, possibly un- 
realistic, efficiency for the thermalization of the released en- 
ergy (Wu, Fabian & Nulsen 2000, Kravtsov & Yepes 2000, 
Bower et al. 2001). 

In this case, alternative sources for ICM heating are 
required, like SN from a primordial star population, the so- 
called Pop III stars (e.g., Loewenstein 2001), and nuclear 
galactic activity associated with QSOs (e.g., Valageas & Silk 

1999, Mc Namara et al. 2000, Yamada & Fujita 2001, Nath 
& Roychowdhury 2002). As for Pop III stars, one expects 
them not to significantly heat the diffuse medium, in order 
not to destroy Ly-a absorbers at lower redshift, z^j 3, and 
not to ovcrpollute the ICM with metals. 

To date, only few attempts have been pursued to ad- 
dress in detail the effects of non gravitational heating on 
X-ray observable quantities with numerical hydrodynami- 
cal simulations (Navarro, Frenk & White 1995, Lewis et al. 

2000, Borgani et al. 2001a, Valdarnini 2002). Bialek, Evrard 
& Mohr (2001) have realized a large set of moderate reso- 
lution simulations of clusters and groups, setting different 
entropy floors at very high redshift. They found that agree- 
ment with observational constraints on X-ray cluster scaling 
properties requires an entropy floor of about 100 keV cm 2 , 
which, at their heating redshift Zh ^ 21, correspond to an 
extra energy of about 0.2 keV/part. 

As an alternative to extra-heating, radiative cooling has 
been also suggested as a mechanism to break ICM self- 
similarity. In this case, X-ray luminosity is suppressed as 
a consequence of the reduced amount of gas in the hot dif- 
fuse phase (Pearce et al. 2000, Muanwong et al. 2001). In 
addition, low-entropy gas in central cluster regions is selec- 
tively removed from the hot phase as a consequence of its 
short cooling time, thus leaving only high entropy gas (Voit 
& Bryan 2001). However, cooling in itself is known to con- 
vert into a cold stellar phase a large fraction of gas, > 30% 
(e.g., Lewis et al. 2000, Dave et al. 2001), much larger than 
the ^ 10% baryon fraction locked into stars indicated by 
observations. This underlines the need for extra heating to 
prevent overcooling of the ICM (e.g. Suginohara & Ostriker 
1998, Prunet & Blanchard 1999, Balogh et al. 2001). 

The main aim of this paper is to study in detail the ef- 
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feet of non-gravitational energy injection into the ICM and 
how much of this energy is needed to correctly reproduce 
different X-ray observable properties of groups and clusters 
of galaxies. In a previous paper (Borgani, Governato, Wad- 
sley et al. 2001a, Paper I hereafter) we concentrated on the 
effect of non-gravitational heating on the entropy pattern of 
the ICM. As a main conclusion, we showed that the entropy 
excess in the central regions of poor clusters (Ponman et al. 
1999) requires an extra heating of about 1 keV/particle. This 
paper presents an extended analysis of a larger set of simu- 
lations, which include several schemes of ICM extra heating. 
We will also show results from a simulation of a poor group, 
which includes gas cooling and a star formation algorithm to 
remove dense cold gas particles from the SPH computation. 
This simulation will be used to discuss the role that cooling 
plays in determining the ICM thermodynamics. 

The plan of the paper is as follows. In Section 2 we 
give a short description of the Tree+SPH code Gasoline and 
provide details about the simulations. Section 3 describes 
the different recipes for non-gravitational heating that we 
used. Section 4 is devoted to the presentation of the results. 
In Section 5 we discuss the effect of cooling, while a summary 
of our main results and our conclusions are presented in 
Section 6. 



2 NUMERICAL METHOD 
2.1 The Code 

We use Gasoline, a new N-body/Smoothed Particle Hydro- 
dynamics code. Gasoline performs gravity and hydrodynam- 
ics operations using a binary tree-structure with periodic 
boundary conditions for cosmology. It uses concurrent mul- 
tiple timesteps for increased throughput and runs in paral- 
lel on a wide variety of architectures using, e.g., pthreads, 
shared memory or the Message Passing Interface (MPI). 
The code handles gravity, pressure gradients, hydrodynamic 
shocks, radiative cooling, photoionizing UV background and 
was extended for this work to incorporate external heat- 
ing sources. For a full description of Gasoline, including ex- 
tended tests of the code, see Wadsley, Stadel & Quinn (2002, 
preparation). For a description of the gravity code on which 
Gasoline is based, see Stadel (2001). The external heating 
SN described in Section 3 is treated within the code as an 
additional component added locally to the standard SPH 
compressive heating ("PdV") term. 

As a demonstration of the effectiveness of Gasoline for 
cluster simulations, we evolve the initial conditions from the 
"Santa Barbara Cluster Comparison" (Frenk et al. 2000). 
We resolve the resulting luminous X-ray cluster with ~ 
15000 gas particles within the virial radius (2.7 Mpc) at the 
current epoch with 64 3 gas and 64 3 dark matter particles 
in the whole box. This resolution was used by the majority 
of the cluster comparison SPH codes including the widely 
used HYDRA code (Couchman 1991). The test evolved the 
gas adiabatically without additional heat sources or sinks. 
The gas profiles for Gasoline with curves for HYDRA and 
the average result from the cluster comparison are shown 
in Figure hi The innermost radial bin plotted in each case 
is indicative of the resolution. The variation seen is simi- 
lar to that found between cluster comparison code results, 



even among different SPH codes, due to differing implemen- 
tations and the interpretation of the initial condition data 
which were given to the coders at 256 3 resolution. The cen- 
tral temperature in particular is sensitive to shock treatment 
and, most of all, minor variations in the timing of small sub- 
clump mergers (for an extended discussion see Wadsley et al. 
2002, in preparation). Except for these small variations, the 
Gasoline results are in excellent agreement with the other 
codes. 

2.2 The Simulations 

We simulated at high resolution four halos taken from a 
DM-only cosmological simulation within a 100 Mpc box of 
a ACDM model with fi m = 0.3, fi A = 0.7, crg=l and Hubble 
constant Ho — 70 km s" 1 Mpc -1 . The baryon density pa- 
rameter was chosen to be fit, = 0.019 h~ 2 , as inferred from 
the deuterium abundance detected in QSO absorption lines 
(e.g., Buries & Tytler 1998), which corresponds to a baryon 
fraction /bar — 0.13. 

From the simulated box we selected four halos with 
masses corresponding to virial temperatures in the range 
0.5< T< 8 keV. With this choice, we are sure to encom- 
pass the temperature interval where the Lx~T relation is 
observed to steepen substantially and where the effect of 
energy injection should be more important. The sizes of the 
four simulated structures are of the same order as those 
of the Coma cluster (e.g. Geller, Diaferio & Kurtz 1999), 
the Virgo cluster (Bingelli, Tammann & Sandage 1987), the 
Fornax group (e.g. Drinkwater et al. 2000) and of a smaller 
Hickson-like group (e.g. Pildis, Bregman & Evrard 1995). In 
the following we indicate with these names the correspond- 
ing simulated structures. 

Using the so-called renormalization technique (Katz & 
White 1983), the simulation within the periodic box is cen- 
tered around each re-simulated halo, the mass distribution 
and the perturbation spectrum are resampled at higher res- 
olution within the Lagrangian region that contains the halo 
of interest. The high-resolution region extends typically for 
a few comoving Mpc around the halo. The outer regions are 
resampled at decreasing resolution. As the four halos span 
almost two decades in mass, the numerical resolutions differ 
accordingly, so as to keep comparable the number of parti- 
cles within the virial radius and the force resolution in units 
of the virial radius. Typically 0.5-1.5 million particles are 
used for each run, with about half of them to simulate the 
high resolution central regions. A summary of the character- 
istics of the simulations are reported in Table |l| Note that 
halo particle number refers to the number of particles within 
the virial radius, and not just in the high resolution region, 
which typically contains a few times more particles. In or- 
der to investigate the effect of resolution, the Virgo, Fornax 
and Hickson runs have been realized with different mass and 
force resolution. The high-resolution version of each run has 
a force softening which is at most 1% of the virial radius, 
while the mass resolution is chosen so that at least 5 x 10 4 
gas particles fall within that radius by z — 0. For all the 
runs we used r\ — 0.2 and a force accuracy parameter of the 
tree algorithm of 8 = 0.7 (9 = 0.5 for z > 2), as suggested 
by Moore et al. (1998) and later verified by Power et al. 
(2002), for the parameters regulating time accuracy of the 
integration and the accuracy of the tree algorithm. 
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Figure 1. "Cluster comparison" gas profiles. The Gasoline results at 64 3 (diamonds) for the temperature (left) and X-ray luminosity 
(right) profiles are in excellent agreement with the averaged results from different codes (solid line) and HYDRA (plus symbols) profiles 
from the cluster comparison paper (Frenk at al. 2000). 



Table 1. Characteristics of the simulated halos and numerical parameters of the simulations. Column 2: simulation version (HR high 
resolution; LR: low resolution); Column 3: total mass within the virial radius at z=0 (10 13 Mq); Column 4: virial radius (Mpc); Column 
5: mass of gas particles (10 8 Mq); Column 6: mass of DM particles within the re-simulated region (10 8 Mq); Column 7 and 8: number of 
gas and DM particles within the virial radius; Column 9: softening parameter (kpc; runs marked with * have been run also with a twice 
as large softening); Column 10: starting redshift. 



Run 






Rvir 




m DM 


Ngas 


AW 


c 




Coma 




133.6 


2.86 


17.7 


119.2 


0.99e5 


0.90e5 


7.5 


19 


Virgo 


HR 


30.4 


1.75 


2.21 


14.9 


1.70e5 


1.79e5 


7.5 


69 




LR 






17.7 


119.2 


0.21e5 


0.22e5 


15 


49 


Fornax 


HR 


5.91 


1.01 


0.65 


4.41 


1.06e5 


1.18e5 


2.5* 


89 




LR 






5.20 


35.3 


1.33c4 


1.48e4 


7.5 


69 


Hickson 


HR 


2.49 


0.76 


0.65 


4.41 


4.93e4 


4.86e4 


2.5* 


89 




LR 






5.20 


35.3 


0.62e4 


0.61e4 


7.5 


69 



We show in Figure g the gas density maps for the high- 
resolution regions of the simulated structures at z = 0. The 
quite high resolution achieved allows us to resolve numerous 
substructures. Some of them retain part of their individual 
gas content, which survives for a few crossing times before 
being stripped, as well as discontinuities associated with bow 
shocks (see (see Paper I for a discussion of these features as 
seen from the entropy maps of simulations) . It is tempting to 
associate such structures to the small-scale features revealed 
by Chandra observations, which in some cases are associated 
with bow shocks or cold fronts from merging structures (e.g., 
Markevitch et al. 2000, Mazzotta et al. 2001, Ettori & Fabian 
2001, and references therein). At z = the Coma-like halo is 
undergoing several merging events with massive sub-groups, 
at different stages of advancement, much like what is seen 
in the XMM observations of the real Coma cluster (Briel et 
al. 2001). 



2.3 Definitions of observables 

The bolometric luminosity for a set of N ga s gas particles is 
defined as 

N ga3 

Lx = i2ssl y -^A c ( Ti ), (i) 

i — l 

where m p is the proton mass, m gas is the mass of a gas parti- 
cle (see Table |l]), pt and Ti are the density and temperature 
at the position of the i-th particle, respectively, and /i = 0.6 
is the mean molecular weight for a primordial gas composi- 
tion with 76% mass provided by hydrogen. Assuming that 
the main contribution to the X-ray emission is from free- 
free bremsstrahlung, the cooling function turns out to be 
A C (T) = 1.2 x 10" 24 (fcT/keV) 1/2 erg cm 3 s" 1 (e.g., Navarro 
et al. 1995). 

Under the assumption of isothermal gas following the 
spherically symmetric DM distribution, it can be shown that 
the bolometric X-ray luminosity is 



Figure 2. The map of the gas density for the HR version of the four simulated structures. The Coma, Virgo, Fornax and Hickson runs 
correspond to the upper left, upper right, lower left and lower right panels, respectively. Each box correspond to a physical size of 10 
Mpc. 
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(e.g., Eke et al. 1998) where f ga s is the gas fraction within 
the cluster, A„; r the overdensity corresponding to virializa- 
tion (~ 178 for Q m = 1; Eke, Cole & Frenk 1996), p the 
average cosmic density and M V i r the virial cluster mass (i.e., 
the mass enclosed within the radius R V i r within which the 
average density is A vir p). Under the above assumptions, 
mass and temperature are related according to 



k B T 



1.38 



Mvir 



1/3 
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(e.g., Eke et al. 1996). The /3 parameter is defined as the 
ratio of the specific kinetic energy of the collisionless DM to 
the specific thermal energy of the gas, 



pm p al 
k B T 



(4) 



being a v the one-dimensional velocity dispersion of DM par- 
ticles. For a Navarro-Frenk- White (1997, NFW) density pro- 
file, F(c) is a function of the concentration parameter c only: 



F(c) = (? 



1 - (1 + c ) 



ln(l + c) - c/(l + c)] 2 



(5) 



For the sake of comparison between the prediction of eq.(g|) 
with results from simulations, we compute the concentration 
parameter as a function of the halo mass for our choice of 
cosmological parameters following the prescription of Eke, 
Navarro & Steinmetz (2001). 



6 Non-gravitational gas heating 



Besides the bremsstrahlung emission, line emission gives 
a non-negligible contribution at temperatures T < 2 keV. 
To account for it and correctly compare our results to ob- 
servations, we use the Raymond-Smith (1977) code in the 
XSPEC package to estimate the emissivity of a plasma with 
metallicity Z — O.SZq. 

We define the mass-weighted temperature of the cluster 

as 



— Y 

N gas ^ 



T, 



(6) 



where the sum is over all the particles falling within the clus- 
ter virial radius. A more useful quantity for comparison with 
observations is the emission-weighted temperature, which is 
defined as 



n3/2 



1/2 



(7) 



This definition strictly holds only for bremsstrahlung emis- 
sivity. We verified that final values of T em are essentially 
unchanged if we account for the contribution from metal 
lines. 

Finally, we define the gas entropy carried by the i-th 
particle as 

m 



n. 



(8) 

where n e> j is the electron number density associated with 
that particle. 

In Table ^| we give the values of bolometric luminosity, 
mass and emission-weighted temperatures for the simulated 
structures and for the different pre-heating schemes (see be- 
low). 



2.4 Effects of resolution 

It is generally regarded that a few thousands of DM parti- 
cles within the virial radius are required to describe global 
quantities, such as total virial mass M V i r and velocity dis- 
persion a v , and as many gas particles are also sufficient to 
correctly estimate the ICM temperature. However, a resolu- 
tion higher by about one order of magnitude is required to 
correctly describe the X-ray luminosity: since this quantity 
depends on the square of the local gas density, its correct 
computation requires the details of cluster structure in the 
innermost regions and gas clumping to be accurately re- 
solved (e.g., Anninos & Norman 1996, Navarro et al. 1995, 
Eke et al. 1998, Bryan & Norman 1998, Lewis et al. 2000). 

Resolution is even more an issue if one wants to properly 
resolve the internal structure of the ICM, including shock 
patterns of the infalling gas, which carry information on the 
physics of the diffuse medium. The higher the resolution, 
the larger the number of substructures surviving within the 
virial radius of the main cluster halo. The largest sub-halos 
are able to keep a fraction of their gas and can carry it down 
to the center of the main halo, possibly shocking the gas 
and raising the central entropy. The X-ray luminosity pro- 
files also depend on the underlying matter profile, for which 
about a hundred thousand particles are necessary to avoid 
underestimating its central density (Moore et al. 1998). 

In Figure H we compare the relation between X-ray 
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Figure 3. The effect of numerical resolution and line— emission 
on the relation between X— ray luminosity and mass-weighted 
temperature for the GH runs. Open and solid circles refer to 
pure bremsstrahlung luminosity for the set of low-resolution and 
high-resolution runs, respectively (see Table 1). Squares are for 
luminosity computed by accounting for the contribution from 
line emission, estimated from a Raymond-Smith code assuming 
Z = 0.3Zq for the global ICM metallicity. For the Fornax and 
Hickson runs, the almost overlapping lower and upper circles are 
for softening parameter e = 5 and 2.5 kpc, respectively (see text). 
The dashed line is the prediction of the model described in the 
text, based on a bremsstrahlung cooling function, an NFW profile 
for the dark matter halo, with concentration parameter appropri- 
ate for the simulated cosmology, and an Af v j r — T relation given by 
eq.(§) with (3 = 1. 



luminosity and mass-weighted temperature, T mw , for the 
HR and LR runs with gravitational heating (GH) only. In- 
creasing the resolution has a non-negligible effect on the 
estimated X-ray luminosity, while, as expected, has only a 
marginal effect on the mass-weighted temperature. The Lx 
value for the HR run of the Virgo cluster is ~ 25% higher 
than for the LR run, the difference increasing to ~ 50% for 
the Fornax group and to ~ 100% for the Hickson group. We 
also verified that decreasing the softening by a factor two 
for the Fornax and Hickson runs, while keeping the mass 
resolution fixed, increases Lx by about 10%, thus showing 
that the adopted spatial resolution is adequate to resolve 
all the structures which are responsible for the X-ray emis- 
sion. In fact, at the highest resolution of our runs, the num- 
ber of massive subhalos (i.e., with circular velocity larger 
than one-tenth of that of the main halo) is likely to have 
converged (Ghigna et al. 2000). We expect spatial resolu- 
tion to be even less of an issue when considering simulations 
with extra heating (see below). Indeed, in this case the gas 
is put on a higher adiabat and, therefore, does not follow 
the small-scale potential wells which are characterized by a 
virial temperature of a few tens of keV. 

We stress that the HR runs predict a Lx~T relation 
which agrees well, both in slope and normalization, with 
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Table 2. Quantifying the simulated ICM at z = 0. Column 1 shows the type of the run and corresponding label. For the runs with 
SN feedback r\ represents the fraction of the total energy released by SN which is thermalized into the ICM. For the S-50 and S-100 
runs, Sfl indicates the value of the entropy floor set at z = 3. Lf, r is the X-ray luminosity from bremsstrahlung emission (10 43 ergs" 1 ); 
Lrs is the X-ray luminosity using a Raymond-Smith code to account for the contribution from line emission, assuming Z = 0.3Zq for 
the average gas metallicity 10 43 ergs _1 ); T mw and T ew are mass— weighted and emission-weighted temperature within R v i r (keV). The 
asterisks mark the simulations for which only the LR runs are available. 
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Figure 4. Effect of resolution on the profiles of gas density (upper panels) and specific entropy (lower panels) for the GH version of the 
Virgo and Hickson simulations. In each panel continuous and dashed curves correspond to the HR and LR runs, respectively. Curves 
become thinner at the radius containing 100 gas particles. The vertical arrows mark the scale corresponding to three times the value of 
the Plummer-equivalent force softening for the LR and HR runs (see Table 1), therefore close to the actual spline softening length. The 
dotted lines in the lower panels show the S oc R 11 entropy profile predicted by the semi-analytical model by TN01. 



the prediction of the scaling model of eq.(g) for a pure 
bremsstrahlung cooling function. We consider this agree- 
ment as a convincing indication that the HR runs provide 
a correct description of the gas distribution, which requires 
using at least ~ 5 x 10 4 particles within the virial radius. 
This result is in agreement with that from resolution studies 
involving the collisionless component only (e.g. Moore et al. 



1998), which established the minimum number of particles 
required to correctly model the central profile of dark matter 
halos. Fig. ^ also highlights a common problem of simula- 
tions of large cosmological volumes at fixed mass resolution: 
as resolution becomes worse for smaller masses, Lx becomes 
underestimated. For instance, taking the same mass and spa- 
tial resolution for the Virgo, Fornax and Hickson runs would 
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produce spurious steepening of the Lx~T relation. Finally, 
while a pure bremsstrahlung emissivity is a good approx- 
imation at T> 2 keV, the effect of line emission is shown 
to become important for smaller/cooler systems. The net 
effect of properly accounting for them is that of flattening 
the Lx~T relation, thus further increasing the discrepancy 
with respect to observations at the scale of poor clusters and 
groups. 

In Figure ^ we show the effect of the resolution on the 
profiles of gas density and mean entropy per particle. There 
is an overall good convergence of these quantities, at the 
radii encompassing at least 100 gas particles (tick curves), 
with slight deviations at somewhat larger scales for the Hick- 
son run. This number of particles corresponds to about three 
times the number of neighbors over which the SPH smooth- 
ing length is estimated. Therefore, the radius containing 
this many particles indicates the smallest scale which is 
resolved by the SPH kernel with the given softening. On 
smaller scales the gas density tends to be somewhat under- 
estimated and so is the X-ray luminosity, while two-body 
heating causes some flattening of the entropy profile (e.g., 
Steinmetz & White 1997, Yoshida et al. 2001). As mass and 
spatial resolutions are increased, gas and entropy profiles 
steepen in the innermost regions. The absence of an entropy 
core agrees with predictions from analytical models based 
on gravitational gas heating and spherical accretion: gas re- 
siding in the central cluster regions has been accreted at the 
very beginning and, therefore, never significantly shocked, 
thus preserving its initial low entropy. 

Besides extra heating, a further possibility to increase 
the entropy level in central cluster regions could be by 
anisotropic accretion of gas which has been previously gravi- 
tationally shocked by large-scale filaments. In order to verify 
this possibility, we identify at z = gas particles in the cen- 
tral region of a halo and trace them back at higher redshift. 
As an example, we show in Figure ^ the density-entropy 
scatter plot for the gas particles falling at z — within 
0.1J?„i r in the HR Virgo simulation, and their correspond- 
ing distribution at z = 1. In the left panel these particles 
define a darker strip, which correspond to the cluster core, 
while the lighter strip at lower entropy witnesses the pres- 
ence of a merged subclump which still preserves its identity 
(see also Fig. El). At z = 1 these same particles define sev- 
eral individual subclumps, which by the present time have 
merged together at the center of the main halo. Such groups 
at z — 1 correspond to lower entropy structures, in line with 
the expectation that the gas within smaller structures un- 
dergoes weaker shocks. Gas possibly accreted along filaments 
would correspond at z = 1 to structures with density con- 
trast 5 ~ 10 containing shocked gas with S^t 50 keV cm 2 . 
Actually, there is no trace of such gas particles in the right 
panel of Fig. |B| This shows that gas accreted along filaments 
does not penetrate efficiently into the cluster central region 
and contribute to increasing the entropy there. 



3 BEYOND GRAVITATIONAL HEATING 
3.1 Entropy floor 

Our first scheme for non-gravitational heating is based on 
setting a minimum entropy value at some pre-collapse red- 
shift (e.g. Navarro et al. 1995, Bialek et al. 2001, TN01). For 
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Figure 5. Gas entropy vs. overdensity for particles of the Fornax 
group which fall within 0.1i? v ir at z = (left panel) and traced 
back to z = 1 (right panel). 



gas with local electron number density n e and temperature 
T, expressed in keV, at redshift z, we define the entropy as 

-2/3 



T 

2/3 



/ba 



l + X 



p(z h )(l + S g 



TkeVcni ,(9) 



where p(z) — po(l + z) 3 is the average cosmic matter density 
at redshift z, S g the gas overdensity and X the hydrogen 
mass fraction. We choose two values for this entropy floor, 
Sfl = 50 and 100 keV cm 2 , that bracket the observed values 
in small groups and clusters (Ponman et al. 1999, Lloyd- 
Davis et al. 2000). 

We assume Zh — 3 for the reference heating redshift, 
since it is close to the epoch at which sources of heating, like 
SN or AGNs, are expected to reach their maximum activity 
(see below). At z = 3, we select all the gas particles with 
overdensity 5 g > 5, so that they correspond to structures 
which have already undergone turnaround. After assuming a 
minimum floor entropy, Sr, each gas particle having Si < S% 
is assigned an extra thermal energy, so as to bring its entropy 
to the floor value. We estimate the amount of energy injected 
in the ICM in these pre-heating schemes by selecting at 
z = all the gas particles within the virial radius and tracing 
them back to z = 3. We find that taking Sn — 50 keV 
cm 2 amounts to give an average extra heating energy of 
Eh = %Th — 1.4 keV/part for particles that end up within 
the virial radius of the Fornax and Hickson groups at z = 0, 
and Eh — 0.9 keV/part for the Virgo cluster and Eh — 0.8 
keV/part for the Coma cluster. Such values are twice as 
large for Sn = 100 keV cm 2 . 

In the semi-analytical model by TN01, the gas is as- 
sumed not to have undergone any significant gravitational 
shock heating before Zh, so that eq.(t}|) actually gives the 
extra energy to be provided to the gas particles to bring 
them to the appropriate adiabat. The hierarchical clustering 
scenario, instead, predicts that a significant amount of non- 
linear structures already exist at Zh = 3, with gas particles 
raised by gravitational shocks to an entropy level already 
higher than the floor to be set. In our heating scheme, the 
thermal energy of such particles is left unchanged. 

This heating scheme can be defined as an external one, 
in the sense that it is also targeted at particles not belong- 
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ing to collapsed structures. In principle this may violate con- 
straints on the entropy of the diffuse high-z IGM, as inferred 
from observations of Ly-a absorption systems. If we identify 
a Ly-a system as a structure with typical gas overdensity 
<5 g ~ 30, it has at z — 3 a temperature T ~ 2.5 x 10 -3 keV, as 
determined by the condition of ionization equilibrium with 
the UV-background (e.g. Haardt & Madau 2001). This cor- 
responds to a typical entropy S — 5-10 keV cm -2 , thus 
about an order of magnitude smaller than that relevant for 
ICM heating. Therefore, an external heating such as the one 
imposed in our simulations would destroy Ly-a absorbers. 
In this sense, setting an entropy floor must just be consid- 
ered as a approach guided by ICM phenomenology rather 
than motivated by general astrophysical arguments. 

Even in this simple heating scenario, a further param- 
eter is represented by the epoch at which the entropy floor 
is created. As the heating redshift is changed, two opposite 
effects compete in determining the energy budget required 
to create a given entropy floor. On the one hand, the higher 
Zh, the smoother the particle distribution, the smaller the 
number of particles at high 8 g , which require a lot of heating 
to increase their entropy to the desired level. On the other 
hand, the higher z^ , the higher the overall cosmic mean den- 
sity p(z), the smaller the amount of heating energy needed 
at fixed 5 9 . We check the effect of changing z^ by also run- 
ning low-resolution simulations of the Fornax group with 
Zh = 1, 2 and 5. We will discuss in the following the effect of 
changing z^ on the final results. The heating energies cor- 
responding to the different zh are shown in the top panel of 
Figure |l3| . 

3.2 Energy feedback from supernovae 

This pre-heating scheme is based on computing the star- 
formation rate within clusters as predicted by a semi- 
analytic model of galaxy formation, a technique first in- 
troduced by White & Frenk (1991), Kauffmann, White & 
Guiderdoni (1993) and Cole et al. (1994), and subsequently 
adopted and developed by several authors (e.g. Somerville & 
Primack 1999, Wu, Fabian & Nulsen 2000, Cole et al. 2000). 
Here we use a variation of the scheme described by Menci 
& Cavaliere (2000), and we refer to this paper for a detailed 
description of the method. 

In our approach, the merging history of dark-matter 
halos, having the same mass of the simulated ones, is pre- 
dicted according to the Extended Press-Schechter theory 
(EPST, e.g. Lacey & Cole 1993). The processes of gas cool- 
ing, star formation and stellar feedback within galaxy-sized 
structures are described by means of a suitable parametriza- 
tion of the corresponding star formation rate and reheated 
mass, which is ejected back into the hot diffuse medium. 
The free parameters of the model are chosen so as to re- 
produce observed properties of the local galaxy population 
(Tully-Fisher relation, B- and K-band luminosity functions, 
disk-sizes) . 

In Figure we show the integrated star formation his- 
tory m»(z, Mo) of all the condensations which are incorpo- 
rated into a structure of total mass Mo by the present time. 
The curves in Fig. ^correspond to Mo values, which span the 
whole range from galaxies to rich clusters. As Mo grows the 
local star formation rate (SFR) decreases, while the value 
at 2 increases and its peak is attained at higher redshift. 
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Figure 6. The integrated star formation rate for the condensa- 
tions ending up in a present-day structure of mass Mo is plotted 
as a function of redshift for the different values of Mo shown in 
the labels. For comparison, the dotted curve shows the cosmic 
star formation rate predicted by our semi-analytical scheme. 



This is expected in the hierarchical clustering picture, since 
more massive structures are originated from more biased 
regions, where initially the star formation rate was higher. 
The large consumption of cold baryons in the progenitors of 
Mo at large z results in a smaller amount of available star- 
forming gas at small z, yielding the sharply declining ro, at 
z->0. 

After assuming a Salpeter initial mass function 
(Salpeter 1955), the SFRs are then used to derive the SN 
heating appropriate to each system that we simulate, based 
on its virial mass. During the evolution of each halo, this 
energy is shared among all the gas particles having S g > 50. 
This overdensity threshold, which roughly corresponds to 
the density contrast at the virial radius, guarantees that 
gas heating takes place inside virialized regions. We verified 
that the final results do not change if only gas particles with 
S g > 500 are heated. Under the extreme assumption that 
all the energy released by SN is thermalized into the ICM 
(i.e. rj = 1 for the SN efficiency) this scheme dumps a total 
amount of about 0.35 keV per gas particle. For the Hickson 
group we also run two more simulations with rj — 0.1 and 
rj — 4. In the first case, the energy budget is so small that 
the final results are essentially indistinguishable from the 
GH run (see Table ^) and, therefore, we will not comment 
further on this low-efficiency case. As for the second case, it 
corresponds to an amount of extra energy similar to that of 
the S-50 entropy floor. This allows us to check whether the 
final results depend on the way a fixed amount of energy 
is dumped into the diffuse medium. Such a large amount 
of energy can be interpreted as associated with any heating 
source whose evolution follows that of the SFR. This is likely 
the case for AGNs, whose emissivity per unit volume has 
been suggested to evolve according to the SFR (e.g. Cava- 
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liere & Padovani 1989, Boyle & Terlevich 1998, Franceschini 
et al. 1999). 

A word of caution should be mentioned here about our 
implementation of SN heating. We assume that all SNe re- 
lease a fixed amount of energy, Esn = 10 51 ergs. In the 
so-called hypernovae scenario, energies larger by at least 
one order of magnitude may be expected (e.g. Loewenstein 
2001). 

Furthermore, in our scenario only Type II SN are in- 
cluded, while Type la SN may in principle give a significant 
contribution to the ICM energetics (e.g., Pipino et al. 2002). 
Also, most of our simulations assume no radiative losses in 
SN explosions, so that the released energy is all thermalized 
in the ICM, which is a reasonable approximation if most 
of the SN explode in a hot and rarefied medium. Finally, 
the feedback used in the semi-analytical modeling of galaxy 
formation is not the same as that adopted for heating gas 
particles in the simulations. For all these reasons, we tend to 
consider our SN heating scheme as a realistic but approxi- 
mate recipe to estimate the amount of energy expected from 
a stellar population, rather than a self-consistent approach 
to including the effect of SN explosions on the ICM energet- 
ics. 



4 RESULTS 

4.1 Gas— density and temperature profiles 

A unique signature of the lack of self-similarity in galaxy 
clusters is provided by the different profiles of gas density 
observed for clusters of different temperatures (e.g. Ponman 
et al. 1999, Lloyd-Davis et al. 2000), with colder systems 
having shallower internal profiles. In the upper panels of Fig- 
ure ^ we show how extra heating can affect the gas profiles. 
In the absence of any extra heating, all the simulated struc- 
tures display the same gas density profile. This result agrees 
with the expectation of the self-similar scaling paradigm and 
confirms that non-linear gravitational processes, like accre- 
tion shocks, do not introduce any characteristic scale. The 
only structure slightly deviating from self-similarity in the 
GH runs is the Coma cluster. Looking at its evolution, it 
turns out that this structure underwent a very recent ma- 
jor merger event, which redistributed the gas in the central 
regions. Therefore, its slightly softer gas profile is just the 
signature of this recent merger. 

Self-similarity is broken as additional heating is in- 
cluded. SN feedback has some effect only on the gas profile 
for the Hickson group, our less massive halo. This is in line 
with the expectation that ICM thermodynamics ought to be 
changed only when the non-gravitational heating is compa- 
rable to the gravitational one. In fact, SN heating provides 
~ 1/3 keV per gas particle, which is similar to the mass- 
weighted temperature of the Hickson group. As for the S-50 
runs, the entropy floor corresponds to a large enough heat- 
ing that only the Coma cluster remains almost unaffected. 
This result is in line with that based on the shape of entropy 
profiles presented in Paper I. 

A relevant quantity which is directly related to the gas 
profile is the baryon fraction, /bar, within a given radius. 
The measurement of the baryon fraction in clusters, in com- 
bination with constraints on f2j from primordial nucleosyn- 
thesis, is considered as one of the fundamental tests for the 



density parameter Q m (e.g., White et al. 1993, Evrard 1997, 
Ettori 2002). We show in Figure ^|the radial dependence of 
/bar for the Virgo and Hickson runs with the different heat- 
ing schemes. If gas is only subject to gravitational heating, 
the cosmic baryons fraction is recovered at about 0.3-0.4 
J?vir, which is comparable to the effective scale where / gas 
is determined from X-ray data. As for the Virgo cluster, 
the S-50 scenario would lead only to a ~ 10-15% underes- 
timate of the correct value at about half virial radius. The 
situation is quite different for the Hickson group, whose gas 
fraction is diminished out to the virial radius by heating 
from the entropy floors. This provides a warning about the 
determination of / gas and the inference of f2 m based on poor 
clusters and groups. The heating from the SN scheme with 
400% efficiency (SN-4) produces a /g aa profile quite similar 
to the S-100 scheme within the virial radius. Quite inter- 
estingly, the S-100 scheme involves an extra-heating energy, 
Eh ~ 2.8 keV, which is almost twice as large than for the for- 
mer scheme. This indicates that a redshift-modulated heat- 
ing is more efficient, in terms of the required energy budget, 
than an impulsive pre-collapse heating. 

As for the three-dimensional temperature profiles (lower 
panels of Fig. [?]), the gas looks almost isothermal, in the ab- 
sence of any extra heating, out to ~ 0.2J? v ir, while a fairly 
steep gradient appears at larger radii. Dips in the T-profiles 
indicate the presence of merging structures, whose lower 
virial temperature brings colder gas in the main body of 
the simulated system, before being stripped by ram pressure 
and thermalized. Such structures are progressively washed 
out in the pre-heated runs, as a consequence of the smaller 
amount of gas that subgroups are able to keep in such cases. 
The extra heating also steepens the temperature profiles in 
the internal cluster regions, with a more pronounced gra- 
dient for smaller systems. An increase of gas temperature 
in the central regions is the consequence of suppressing gas 
density while keeping unchanged the external gas pressure. 

In the last few years, observational data have reached 
good enough precision to allow for spatially-resolved spec- 
troscopic observations of the ICM for fair number of galaxy 
clusters. Yet, a general agreement about the temperature 
profiles in the external cluster regions, not affected by cool- 
ing, has still to be reached (e.g., Markevitch et al. 1998, 
White 2000, Irwin & Bregman 2001). Based on Beppo-SAX 
observations, De Grandi & Molendi (2002) analyzed pro- 
jected temperature profiles for a set of 21 clusters. They 
found the gas to be isothermal out to ~ 0.2_Ri8o, with a 
fairly steep negative gradient at larger radii. In Figure ^ we 
compare result from our simulations to the average profile 
by De Grandi & Molendi (2002). Since their clusters are all 
quite hot (T> 4 keV), we use only results from the Coma 
run. The simulation results are obtained by averaging over 
the temperature profiles projected along three orthogonal 
directions. The simulation produces a fairly steep gradient 
at R^, 0.2i?i8o, quite similar to the observed one. Forthcom- 
ing data from Chandra and Newton-XMM (e.g. Arnaud et 
al. 2001, Allen, Schmidt & Fabian 2001) will certainly refine 
the determination of temperature profiles in central cluster 
regions, where cooling is also expected to play a significant 
role. There is no doubt that such observational advances 
will challenge the ability of numerical simulations to treat 
the relevant physical processes. 
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Figure 7. Profiles of gas density in units of the virial density (upper panels), and temperature in units of T ew (lower panels) at z = 0. 
Left, central and right panels are for the GH, SN and S-50 runs, respectively. Dotted, solid, short— dashed and long-dashed curves refer 
to the Coma, Virgo, Fornax and Hickson runs, respectively. 
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the SN heating scheme with efficiency raised to 400% (SN-4). The horizontal dashed line indicates the cosmic baryon fraction. 



4.2 The mass—temperature relation 

The M-T relation represents a key ingredient when com- 
paring observational data on the cluster X-ray temperature 
function with the mass function predicted by cosmological 
models. Quite recently, the different results obtained by dif- 
ferent authors on the power-spectrum normalization, as, 
from cluster XTF and XLF have lead to a discussion on 



the correct M-T relation to be used (e.g. Pierpaoli, Scott 
& White 2001, Borgani et al. 2001b, Seljak 2002, Viana, 
Nichol & Liddle 2002, Ikebe et al. 2002). In past years, dif- 
ferent analyses have shown that eq. Q provides a good fit 
to results of simulations including only gravitational heat- 
ing, with values of (3 ranging within the interval 0.9i$ (3^ 1.3 
(e.g. Evrard, Metzler & Navarro 1996, E96 hereafter, Bryan 
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Figure 9. The projected temperature profile for the GH and S- 
50 Coma runs (solid and dashed lines, respectively), compared 
with the average temperature profile for the set of 21 clusters ob- 
served with Beppo-SAX and analyzed by De Grandi & Molendi 
(2002). Simulation results are obtained by averaging the profiles 
obtained by projecting along three orthogonal directions. The 
cluster— centric radius is expressed in units of Riso, which is de- 
fined as the radius encompassing the density 180p cr it . 



& Norman 1998, Frenk et al. 1998). In particular, E96 
pointed out that temperature measurements provide a quite 
accurate determination of the mass at radii where the mean 
cluster density is 500-2500 p C rit • 

Besides investigating the effect of non-gravitational 
heating on the M-T relation, the excellent resolution of our 
simulations allows to better resolve the structure of cluster 
cores, which is relevant for determinations of the emission- 
weighted temperature. Results on the M v i r -T m w relation 
from our simulations (left panel of Figure [To]) closely fol- 
low the relation of hydrostatic equilibrium of eq.(^) with 
a very small scatter and almost no dependence on non- 
gravitational heating. Although extra heating adds internal 
energy to the gas particles, they quickly get back in hydro- 
static equilibrium, with their temperature determined by the 
depth of the gravitational potential well which is established 
by the dynamically dominant DM component. 

As for observational data, the M-T relation has been 
measured in the last few years by different groups us- 
ing ASCA (Horner et al. 1999, Nevalainen et al. 2000, 
Finoguenov et al. 2001, F01 hereafter) and Beppo-SAX data 
(Ettori, De Grandi & Molendi 2002). Such analyses consis- 
tently find that: (a) M oc T 3/2 for T> 4 keV, but with a 
normalization significantly lower than that found by E96 
from simulations (cf. Ettori et al. 2002); (b) a steeper slope 
for colder systems, possibly interpreted as an effect of pre- 
heating. We compare in the right panel of Figure [l^ data on 
the Msoo-Tew relation by F01, which include data on sys- 
tems down to T ew ^, 1 keV, to results from our simulations 
and to the relation found by E96. The somewhat lower nor- 
malization with respect to the E96 results on the group scale 



is likely to be due to our improved resolution. However the 
effect of extra heating is at most marginal and not sufficient 
to reconcile simulations to data, thus indicating that some 
other physical process should be at work in establishing the 
M-T scaling. 

Finoguenov et al. (2001, F01 hereafter) suggest that 
the difference between observed and simulated M-T rela- 
tion is due to the combined effect of pre-heating and the 
effect of formation redshift on the temperature of the sys- 
tem. However, our simulations show that pre-heating has a 
minor impact on this relation. As for the effect of forma- 
tion redshift, it may introduce a bias in the definition of 
the observational data set: observations could tend to se- 
lect fairly relaxed systems, which formed at higher redshift 
and, therefore, are characterized by a somewhat higher tem- 
perature at a fixed mass (e.g. Kitayama & Suto 1996, Voit 
& Donahue 1998). Our simulations define an M-T relation 
with a small scatter, thus suggesting that differences in the 
formation epoch or differences in the current dynamical sta- 
tus among systems should have a small effect. A larger set of 
simulated clusters would be required to properly address this 
point. Ettori et al. (2002) detect a segregation in the M-T 
relation for cooling-flow and non cooling-flow clusters, the 
latter being characterized by a larger scatter. 

TN01 showed that their ICM model, which incorpo- 
rates the effects of pre-heating and cooling, reproduces the 
observed M-T relation. They also find that the predicted 
relation is weakly sensitive to the value of the pre-collapse 
entropy floor. This suggests that cooling should be respon- 
sible for the lower normalization of the relation, through 
the steepening of the temperature profiles in cluster central 
regions. The effect of cooling on the M-T relation will be 
further discussed in Section 5 below. 

4.3 The luminosity— temperature relation 

The observed relation between bolometric luminosity and 
temperature is considered a standard argument against the 
self-similar behavior of the ICM. Bremsstrahlung emissivity 
predicts Lx oc Mpg as T 1,/2 . Therefore, as long as clusters of 
different mass are scaled versions of each other, then the M- 
T scaling from hydrostatic equilibrium gives Lx oc T X {1 + 
z) 3/2 or, equivalents L x oc M 4/3 (l + z) 7/2 for fi m = 1 
(Kaiser 1986, see Eke et al. 1998, for an extension to low- 
Q m cosmologies). As we also discussed in the introduction, 
this prediction is at variance with respect to observational 
evidence of a steeper relation, Lx oc T~ 3 for TjJ 2 keV 
and, possibly, even steeper for colder systems. This result is 
also in line with the observed slope of the Lx~M relation, 
L x oc M a with a ~ 1.8 ± 0.1 (Reiprich & Bohringer 2002). 

The first determinations of the Lx~T relation for clus- 
ters showed that it has a quite large scatter (e.g. David et 
al. 1993, White, Jones & Forman 1997). A significant part 
of this has been recognized to be the effect of cooling: cen- 
tral spikes associated with cooling regions provide a large 
fraction of total X-ray luminosity, so that differences in the 
cooling structure among clusters of similar temperature in- 
duce a spread in the corresponding Lx values. After correct- 
ing for this effect, different authors (Allen & Fabian 1998, 
Markevitch 1998, Arnaud & Evrard 1999) were able to cali- 
brate a much tighter Lx~T relation. Allen & Fabian (1998) 
also noticed that the Lx~T relation for the hottest clus- 
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Figure 11. The relation between X-ray luminosity and 
emission-weighted temperature at z = 0. Results from simu- 
lations are compared to observational data for nearby clusters 
(data points with errorbars) from different authors. The asterisk 
indicates the SN-4 Hickson run, which corresponds to E h = 1.4 
keV/part for the heating energy. 



of temperatures. Heating with the SN-1 recipe provides some 
steepening of the Lx~T relation at the group scale, although 
not enough to reach agreement with data points. A larger 
suppression of Lx is achieved with the S-50 and S-100 heat- 
ing schemes and, for the Hickson group, with the SN-4 run. 
Quite interestingly, the SN-4 heating scheme requires the 
same heating energy, Eh — 1.4 keV/part as the S-50 en- 
tropy floor, but provides a significantly smaller luminosity, 
as a consequence of the lower gas density in the central re- 
gion of the Hickson group (see also Fig. H). Therefore, a bet- 
ter efficiency is reached in this case by gradually dumping 
energy within the virialized regions of the ICM, rather than 
imposing an impulsive pre-collapse heating on the whole 
turn-around region. 

In a similar way, Figure |l^ shows a similar compar- 
ison with observations for distant, 0.5^ 1.3, clusters. 
Independent analyses have confirmed that data on Lx- 
T relation of distant clusters are consistent with a lack of 
evolution (e.g. Mushotzky & Scharf 1997, Donahue et al. 
1999, Delia Ceca et al. 2000, Borgani et al. 2001b, Stanford 
et al. 2002, Holden et al. 2002). Since no data are avail- 
able on T< 1 keV groups in the above z-range, we do not 
include in this comparison results for the Hickson runs. Al- 
though temperature determinations for distant clusters are 
prone to larger errorbars, the results are quite in line with 
what shown in Fig. [ll| the Lx~T relation of distant clusters 
requires pre-heating of the ICM with EfJZ> 1 keV/part. 



ters tends to flatten to the self-similar scaling prediction, 
possibly suggesting that the effect of extra heating becomes 
negligible for such systems. 

Figure [ll] shows the effect of extra heating on the Lx- 
T relation of nearby clusters. The GH runs are clearly at 
variance with respect to data over the whole sampled range 



4.4 The effect of changing the epoch of heating 

Imposing an entropy floor at some pre-collapse redshift is 
a useful approximation, in that it allows the characteriza- 
tion of the ICM evolution by a single quantity which is 
conserved in adiabatic processes. As already mentioned, the 
heating energy required to establish a given entropy floor is 
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Figure 12. The Lx-T relation for distant z 0.5, clusters. Tri- 
angles connected by lines are simulation results for Coma, Virgo 
and Fornax runs at z = 0.5 (filled triangles with dotted lines) 
and 2 = 1 (open triangles with dashed lines). The three groups 
of simulation results at decreasing Lx refer to the GH, S-50 and 
S-100 runs (for reasons of clarity we do not include the case of 
SN feedback). Data points with errorbars refer to z > 0.5 clusters 
from Beppo-SAX and ASCA data (Delia Ceca et al. 2000, open 
circles), Chandra data extending out to z = 1.27 (Borgani et al. 
2001b, Stanford et al. 2002, filled circles) and from the XMM ob- 
servation of the 2 = 1.26 cluster of the Lockman hole (Hashimoto 
et al. 2002, square). 



a non-trivial function of the heating redshift, Zh, through 
the evolution of p and <5 gas (see eq.g). In general, one may 
also expect the resulting global ICM properties to depend on 
Zh - In their simulations, Bialek et al. (2001) heated at a very 
high redshift, Zh ~ 20, and found results on the required Sa 
to reproduce the observed Lx~T relation which are gener- 
ally consistent with ours. This may suggest that the choice 
for zh is not relevant. However, a close comparison between 
their runs and ours is not straightforward: mass and dynam- 
ical resolutions are higher in our runs and, as discussed in 
Section 2, this is likely to affect observable quantities, such 
as Lx, which are sensitive to the details of gas dumpiness. 

As shown in Figure |l3|, the required Eh for the S-50 
Fornax run changes by about a factor 2, zu = 3 being the 
less efficient choice in terms of energy budget. Both T em 
and central gas entropy at 0.1 R v i r are quite stable, with 
variations of at most ~ 10% and ~ 20% respectively. As 
expected, the X-ray luminosity is the most sensitive quan- 
tity to changes of Zh\ it decreases by about a factor two as 
the heating epoch is pushed back in time from Zh = 1 to 
5. Such relatively modest variations suggest that different 
choices for zu induces only moderate changes in the value of 
Sfi required to reproduce the observed Lx~T relation. 



Figure 13. The effect of changing the heating redshift, 2^, for 
the low— resolution S-50 Fornax run. Results are shown for the 
specific heating energy, the emission-weighted temperature, the 
X-ray luminosity and the entropy at 0.1 R v i r (from upper to lower 
panels). 



5 THE EFFECT OF COOLING AND STAR 
FORMATION 

The simulations discussed so far have been performed with 
the aim of understanding in detail the effect of non- 
gravitational heating on ICM observable quantities. How- 
ever, excluding cooling in our simulations is clearly a sim- 
plification, and its inclusion is likely to change some results 
significantly. Cooling times can be significantly shorter than 
the Hubble time in central regions of clusters, thus causing 
a significant fraction of the gas there to cool down and drop 
out of the hot diffuse, X-ray emitting phase (see Fabian 
1994, for a review on cooling processes in clusters) . 

Cooling has been suggested as an alternative to non- 
gravitational heating for breaking ICM self-similarity. As 
gas undergoes cooling in central cluster regions, X-ray lumi- 
nosity is suppressed as a consequence of the reduced amount 
of diffuse hot gas (Pearce et al. 2000), thus possibly pro- 
ducing a steepening of the Lx~T relation (Muanwong et 
al. 2001, Wu & Xue 2002). The decreased pressure support 
causes the recently shocked, higher entropy external gas to 
flow in, thus causing a net increase of entropy for the hot gas 
left in central cluster regions (Bryan 2000). Recently, Voit 

6 Bryan (2001) have argued that only gas with entropy in 
excess of 100-200 keV cm 2 has long enough cooling time to 
remain in the diffuse phase. Accordingly, the entropy excess 
in central regions of poor groups should be interpreted as an 
effect of cooling, rather than of extra heating. 

In order to check this effect, we run a high-resolution 
simulation of the Hickson group, with the SN-1 heating 
scheme, but also including gas cooling. The cooling function 
adopted assumes gas of primordial composition, with zero 
metallicity, and follows the recipe of Wadsley & Bond (1997). 
We also include heating by UV background by adopting the 
latest version of the model by Haardt & Madau (2001). We 
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Figure 14. The map of the gas density for the SN-1 Hickson run 
including cooling and star formation. The box size is the same as 
in Fig. 2. 
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Figure 15. The entropy profiles for the Hickson simulation with- 
out cooling (for the GH and SN-1 cases) and including the effect 
of cooling and star formation in the SN-1 run. In the computation 
of the entropy for the cooling run, only gas particles belonging to 
the hot diffuse phase have been considered (see text). 



follow the recipe by Katz, Hernquist & Weinberg (1992) to 
convert dense cold gas particles into collisionless stars. Fur- 
thermore, when computing properties of the diffuse ICM, we 
exclude gas particles with temperatures T < 3 x 10 4 K and 
overdensities S > 500 (see also Kay et al. 2002). These par- 
ticles are assumed to belong to the cold medium, although 
not yet removed from the gas phase by the star formation 
(SF) algorithm. The effect of including cooling and SF is 
that of collecting gas into dense knots, while suppressing 



ga s density in the diffuse filamentary structures (see Figure 
[14). The resulting density profile in central regions turns out 
to be steeper, oc R~ a , with a ~ 2 than that for a CDM only 
run (see also Treu & Koopmans 2002). 

In Figure [l^ we show the resulting entropy profile, as 
compared to the GH and SN-1 runs, with no cooling. We 
detect some increase of gas entropy inside R V i r , due to the 
infall of shocked external gas. However, the inclusion of cool- 
ing and SF do not produce any entropy plateau at ~ 100 
keV cm 2 . Although the entropy increase goes in the right 
direction, it may be marginal in accounting for the observed 
entropy excess in central group regions (e.g. Ponman et al. 
1999; see also Paper I). A higher entropy could be possi- 
bly attained by assuming a non-negligible metallicity in the 
cooling function, which would cause an increasing efficiency 
of gas cooling. 

Besides suppressing the amount of hot X-ray emitting 
gas, cooling also causes the formation of high-density con- 
centration of star particles in the innermost cluster region. 
Such concentrations deepens the potential well and causes 
temperature there to significantly increase. Therefore, while 
the A" -ray luminosity decreases, the emission-weighted tem- 
perature increases. The effect of cooling in the M-T ew and 
Lx~T relations of the Hickson run is shown with the filled 
square in Figs. 9-11. Such changes actually go in the right 
direction of improving Lx~T relation (see Figure [ll]) and 
possibly solve the discrepancy with the observed Msoo-Tem 
relation (see also Thomas et al. 2002; Voit et al. 2002). 

However, cooling in itself is a runaway process, due to 
the quick increase of cooling efficiency with local gas den- 
sity. As a consequence, if not counteracted by some feedback 
process, a too large fraction of gas can cool down to a colli- 
sionless phase (e.g., Prunet & Blanchard 1999, Lewis et al. 
2000, Dave et al. 2001, Balogh et al. 2001). We see in our 
simulation that while the overall baryonic fraction inside 
R V ir at 2 = coincide with the cosmic value, 37% of the 
gas is locked in the stellar phase. This large fraction is far in 
excess of the ^ 10% fraction indicated by observations (e.g., 
Balogh et al. 2001). Therefore, our SN feedback scheme falls 
short in preventing the cooling runaway. 

Addressing in detail the issue of preventing overcooling 
through a suitable feedback mechanism is beyond the scope 
of this paper. Due to the biased nature of halo formation in 
hierarchical models, the first massive halos to collapse into 
a protocluster region will eventually merge to form its core. 
These massive halos host the progenitors of giant ellipticals 
in today clusters (Governato et al. 2001, Springel et al. 2001) 
and will likely contribute to most of the energy released in 
the ICM. Whatever the astrophysical source for this energy 
is, it should act in a self-regulated way: cooling switches 
on the feedback mechanism and the subsequent energy re- 
lease inhibits too much gas to leave the hot phase. Besides 
SN, also AGN activity has been suggested as a possible way 
of preventing overcooling, with energy released with a duty 
cycle of a few 10 7 years (e.g. Binney & Tabor 1995, Ciotti 
& Ostriker 2001, Quilis et al. 2001, Bohringer et al. 2002). 
If AGN are responsible for the required feedback, then one 
should expect to see some signature of their recent activity 
in the central regions of those clusters whose cooling has 
been recently switched off. A succesfull and self-consistent 
inclusion of cooling and feedback in numerical simulations 
has still to be realized, the main difficulty lying in the treat- 
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ment of physical processes which involve a large dynamical 
range. The physics describing AGN activity or SN explosions 
involve small scales that are out of the reach of current and 
foreseeable simulations of galaxy clusters, thus leaving as 
the only possibility the inclusion of such process as external 
recipes. 



6 CONCLUSIONS 

We presented results from high resolution Tree+SPH sim- 
ulations of galaxy clusters and groups aimed at studying 
the effect of non-gravitational heating of the intra-cluster 
medium (1CM) on observable X-ray properties of galaxy 
systems. Four simulated structures have been chosen so as 
to encompass the temperature range from a poor Hickson- 
like group to a rich Coma-like cluster. Two different kinds of 
non-gravitational gas heating have been chosen: (i) impul- 
sive heating to set an entropy floor at some redshift (zu = 3 
as a reference epoch); (ii) gradual heating, which follows the 
star-formation rate inside the cluster region, as predicted 
by semi-analytical modelling of galaxy formation, so as to 
mimic the effect of SN energy feedback. Our main results 
can be summarized as follows. 

(a) In order for simulations to reproduce analytical pre- 
dictions for the Lx~T relation of gas sitting in hydrostatic 
equilibrium within a NFW potential well, the structure of 
the gas distribution needs to be resolved with at least 5- 
10 x 10 4 gas particles within the virial radius, with a spatial 
resolution of the order of 1-2% of R v i r - 

(b) The observed slope and amplitude of the Lf, i~T ew re- 
lation are well reproduced by setting an entropy floor of 
50-100 keV cm 2 . Assuming negligible radiative losses, SN 
feedback provides ~ 0.35 keV/part and, as such, has a too 
small effect on the Lx~T relation. A four times larger en- 
ergy, possibly consistent with that obtainable from AGN, is 
required to reproduce the Lx~T relation at the scale of a 
Hickson group. 

(c) In the absence of cooling, the M$oo-T ew relation is 
found to be higher by ~ 40% than that observed, over the 
whole sampled T-range, almost independent of the presence 
and amount of non-gravitational heating. Therefore, non- 
gravitational heating in itself is not able to account at the 
same time for the Lb i-T em and the Msoo-Te™ relations. 

(d) The effect of including cooling in the simulation of a 
Hickson group with SN heating is to suppress Lx by about 
40% and increase T em by about 30%. Therefore, the effect 
of cooling goes in the right direction of decreasing the am- 
plitude of the Lx~T and M$oo-T ew relations at the scale of 
groups. However, the runaway nature of cooling is not ef- 
ficiently counteracted by our SN feedback scheme, so that 
37% of the gas drops out of the hot X-ray emitting phase. 
Such a large amount of cold gas is in contrast with observa- 
tional data on the stellar content of galaxy systems. 

The results presented in this paper demonstrate that 
no unique simple recipe exists to explain X-ray observable 
properties of galaxy systems. Supernova heating provides 
a sizeable, but insufficient, fraction of the non-gravitational 
heating energy required. Further feedback energy from AGN 
may be required. Yet, the mechanisms to thermalize the en- 



ergy associated with QSO activity in the diffuse medium has 
to be understood in detail. 

While much further work is needed, our study suggests 
that a combination of different heating sources (SN and 
AGN) and cooling will hopefully be able to reproduce both 
the correct Lx~T em and Msoo-Tem relations, and to avoide 
gas overcooling. A self-consistent treatment of such complex 
physics is mandatory if the numerical description of the ICM 
has to keep pace with the progress in the observational de- 
scription of the hot baryons in clusters. 
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